Topology induced Kondo effect in hydrogenated Pt nanocontacts 
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It is shown that recent experimental data on electronic transport through Pt nanocontacts in the 
presence of hydrogen admit an explanation in terms of topological and electron-electron correlation 
grounds. A model Hamiltonian, which incorporates two orbitals on Pt atoms, a single orbital on 
hydrogens, and the on-site Coulomb repulsion on the H atoms, is solved exactly, and connected to 
Pt leads described by a Bethe lattice. When two weakly coupled H atoms are placed between the Pt 
electrodes transversally to the transport direction (as recently suggested) a Kondo effect related to 
the symmetry of the Pt-H couplings, stabilizes the conductance around one quantum with a single 
channel contributing to the current, in agreement with the experiments. 
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A recent experimental work reports that electronic 
transport through Pt nanocontacts is strongly modified 
in the presence of hydrogen [l|, |^ . Specifically, while in 
Pt nanocontacts the main conductance peak lies close to 
1.5 conductance quanta Qq = 2e^//i, in a hydrogen at- 
mosphere the main peak shows up near one quantum. 
Although it was suggested [l|, y that the latter result 
was due to a H molecule being placed along the nanocon- 
tact axis bridging the last two Pt atoms of the left and 
right electrodes, an alternative explanation has been put 
forward recently 3]. Extensive quantum chemistry cal- 
culations led the authors to conclude that a geometry in 
which two weakly coupled hydrogen atoms were placed 
transversally to the nanocontact axis (see Figljl, was 
more probable. Most of the results described in 1] found 
a reasonable explanation j^], including vibrational fre- 
quencies and the isotope effect, but the observed unit 
conductance to which a single channel apparently con- 
tributes. The aim of this work is to study the strong cor- 
relation effects that one may expect due to the presence 
of hydrogen, the role played by topology, recently empha- 
sized in an analysis |j| of transport through two quantum 
dots in a similar geometry (dots located transversally to 
the nanostructure axis), and the importance of both in 
understanding the experimental results [ij. 

The model Hamiltonian used to describe the Pt-H-H- 
Pt cluster shown in Fig. ^is: 
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The Pt Hamiltonian includes two isotropic (s-like) or- 
bitals, the first representing the 6s Pt orbital, and the 




FIG. 1: Sketch of the geometry of the Pt2ll2 cluster con- 
nected to bulk Pt (described by means of a Bethe lattice, see 
text). The z axis is taken perpendicular to the cluster plane. 



second some average of the five 5d Pt orbitals 

Hpt = ^ Uaal^aia ^' ebhl^hia ^ 

i=l,2;CT 



(2) 



<ij>;cr 



As in the geometry of Fign]the H atoms are rather far ap- 
part |3 we do not include the interaction between their 
respective s-orbitals. On the other hand, the Hamilto- 
nian incorporates the on-site Coulomb repulsion through 
a Hubbard term, 

Hh = ^ ehhl^h^a + Uhhl^hi^hl^h^i (3) 

i=l,2;CT 

Finally the Pt-H interactions are written as, 

Hpt-H ^ J2 i^hahtaj^ + thbhj^bj^) (4) 



TABLE I: Energies of the Pt and H orbitals and the interac- 
tions amongst them (all in eV), as defined in Eqs. (1-4), used 
to describe the Pt2 and the Pt2H2 clusters. Parameters for 
bulk Pt were taken equal to those of the Pt2 cluster. 
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When the cluster is connected to electrodes, the trans- 
mission across the system is given by p : 



T(E) = —Tilth]. 
h 



(5) 



and the conductance is Q = T{Ef), where Ep is the 
Fermi level that will be hereafter taken as the zero of 
energy. In this expression, Tr denotes the trace over all 
orbitals in the cluster and the matrix t is given by 



i = 4/2G(+)ri/2= r^/^G^-^r^/' 



where the matrix TR{h) is given by z(S^ ' 
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being the self-energies of the right (R) and left 
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(L) leads. The Green function is written as. 
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where Gg is the Green function associated to the iso- 
lated Pt2H2 cluster, which is obtained by means of ex- 
act diagonalization using the Lanczos method [3, ISj. 
The bulk of the electrodes is described by means of a 
Bethe lattice with coordination twelve, while the two Pt 
in the cluster are connected to c bulk atoms. As dis- 
cussed below, this effective coordination c is taken as a 
fitting parameter. In order to single out the contribu- 
tion of individual channels we diagonalize the matrix vt. 
This method, hereafter denoted as the Embedded Clus- 
ter Method (ECM) , is exact only as far as the calculation 
of the Green function of the isolated cluster is concerned, 
and does not account for electron-electron interaction ef- 
fects that may extend beyond its bounds. In the Kondo 
regime, many-body effects extend over a length given by 
the Kondo cloud that is proportional to the inverse of the 
atomic Kondo temperature Tk- In this case, as hydrogen 
and Pt interact strongly, we expect this temperature to 
be large enough to guarantee a rather localized cloud. 
Aiming to highlight the role played by the Kondo effect 
we also solve this system using Unrestricted Hartree Fock 
(UHF) that gives an appropriate description for temper- 
atures greater than the Kondo temperature. 

The choice of the model parameters (given in Table 
I) was guided by ab-initio calculations of the Pt2 and 
Pt2H2 clusters. As that calculation indicates that the 
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FIG. 2: Upper panel: transmission as a function of energy 
for the Pt2 cluster. Total transmission (thick continuous line) 
and the contributions of the d i(continuous line) and s (broken 
line) eigenchannels are shown. Lower panel: same for the 
density of states (DOS). 



dxy orbital is the one that more strongly interacts with 
the H orbital, orbital b in Eq. (2) was assumed to have 
that symmetry in the Pt2H2 cluster. The parameters 
of Table I give a bulk Bethe-lattice Density of States 
(DOS) with a d band more than half-filled and a s band 
less than half-filled. All Bethe-lattice and nanocontact 
calculations were carried out by fixing the charge to one 
(two) electron per H (Pt) atom. Parameter c in Eq. ^ 
was chosen to give the experimental conductance in Pt 
nanocontacts with no hydrogen 0, Q ■ 

Fig|21 shows our results for the transmission and the 
DOS in this nanocontact (with c=5.5), which are in line 
with those reported in [3j. Both channels equally con- 
tribute to give a conductance around I.GQq. Due to the 
symmetry of the dxy orbital there is no indirect (hydro- 
gen mediated) s — d coupling between Pt atoms. In ad- 
dition, direct coupling between s and d orbitals on dif- 
ferent atoms is small enough to be ignored in our sim- 
plified model. As a result, s and d orbitals act indepen- 
dently even in the presence of hydrogen. Our results for 
the transmission in the hydrogenated Pt nanocontact are 
depicted in Fig|31 Comparing these results with those 
without H (see Fig|21) wc note that now the transmis- 
sion reaches a maximum (minimum) of around Qq (0) for 
the d (s) eigenchannel. This result is not reproduced by 
the UHF approximation, as shown in the lower panel of 
FigEl The latter resembles the transmission for the Pt 
nanocontact (see Fig|2Jl, but reduced by the effect of scat- 
tering with H orbitals. The local DOS at the H atoms 
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FIG. 3: Left panels: s and d transmissions versus energy in 
the Pt2H2 cluster as calculated with ECM (upper) and UHF 
(lower). Right panels: same but replacing the Pt d orbital by 
an s-like isotropic orbital. 



is shown in the upper panel of Fig01 The ECM result 
shows the Kondo signature at the Fermi level, which is 
absent in the UHF calculation, that describes the solution 
for T > Tk. The comparison of UHF and ECM results 
clearly suggest the presence of strong many-body effects. 
As there is no direct or indirect hybrization between s 
and d bands, the pair of H can be in the Kondo regime 
with each Pt band simultaneously and independently, as 
described in [^ for a one band model. In this regime, 
both hydrogens become strongly correlated (in a state 
of total spin Sh2 = l)i forming a Correlated Impurity 
of Spin 1 (CISl), side connected to the respective chan- 
nel in our model. Then, the conductance in each channel 
has two contributions coming from: i) direct coupling be- 
tween d or s orbitals, and, ii) coupling mediated by the 
CISl. Charge in the the CISl induce a phase shift of tt 
respect to the direct contribution, as can be argued from 
the Friedel sum rule |10l] . and was also verified numer- 
ically by us. Contributions from both ways are similar 
and therefore a destructive interference take place in the 
s channel. In the case of the d channel there is an ad- 
ditional 7r/3 phase shift, because the H are connected 
with different sign to each Pt. This prevents destructive 
interference in this channel. Moreover, in this case the 
dominant contribution to the conductance comes from 
the impurity channel (see Table I) that is in the Kondo 
regime given conductance very close to 1. Results for the 
transmission and local DOS for the model with two s-like 
orbitals (si and S2) in place of s and d orbitals are also 
shown in FigsOl and ^ In this case there is no Kondo 



FIG. 4: Upper panel: Local DOS at the H in the Pi2-H"2 
nanocontact, as calculated by means of UHF and ECM, with 
either two s Pt orbitals or two orbitals of s and d:^y symmetry 
(see text). Lower panel: Conductance and spin correlations 
calculated with ECM as function of the degree of hydroge- 
nization (see text). 



resonance at the Fermi level, as si and S2 orbitals on dif- 
ferent Pt atoms are now coupled through the hydrogens. 
As shown in FigQ] this coupling splits the Kondo reso- 
nance into two peaks above and below the Fermi level. 
As a consequence, only the direct coupling between Pt 
atoms contribute to the conductance and, as interference 
effects do not show up, the two channels give similar con- 
tributions to the conductance (see FigOJ. This is also 
the case of a H2 molecule placed longitudinally to the 
nanocontact axis, as, in this case, each H atom is con- 
nected only either to the left or right electrodes, giving 
an indirect coupling between s and d bands. 

In order to investigate how the conductance changes 
from a situation with no hydrogen up to that where the 
two H atoms are located as in Fig^ we have carried out 
calculations varying the H-Pt parameters as tf^^ = xnthi, 
where i — a,b and xh varying from (H free nanocon- 
tacts) up to 1.5. The most appealing feature of the re- 
sults shown in the lower panel of FigQlis the wide range of 
Xh over which the ECM gives a conductance close to one 
[l2| . This is due to the electron-electron interactions that 
are playing an essential role in this system. The latter 
statement is corroborated by the spin-spin correlations 
also shown in Fig^ When the hydrogen atoms are not 
coupled (or just weakly coupled) to Pt, H-H correlations 
are antiferromagnetic and the total spin of the molecule is 
zero, as occurs in the free hydrogen molecule. In this limit 
Pt-H spin-spin correlation is zero, as it should. However, 
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FIG. 5; Normalized conductance histogram and normalized 
conductance fluctuations for the hydrogenated nanocontact 
versus the conductance through it. The results were obtained 
by using either ECM (upper panel) or UHF (lower panel). 



just when the conductance stabilizes around one, H-H 
spin correlations become strongly ferromagnetic, and the 
total spin of the molecule raises up to 0.8, while the Pt-H 
bond shows antiferromagnetic correlations. These results 
illustrate the many features that this many-body effect 
shares with the standard Kondo effect of spin 1 already 
obtained in a two parallel quantum dot system |j]. 

The method used in Ij to estimate the number of chan- 
nels that contribute to the current consists of calculating 
the conductance fluctuations obtained in a large num- 
ber of experiments. The experimental results indicate 
that the derivative of the conductance G with respect 
to the bias voltage dG/dV varied randomly obeying a 
bell-shaped distribution whose standard deviation was 
calculated in Q, 



'^GV oc Y^ r^(i-T„) 

V n=l,..,N 



(8) 



where r„ is the transmission probability of channel n and 
N is the total number of channels. We have calculated 
the conductance fluctuations by randomly varying xh 
over the range [0,1.5]. Couplings between Pt atoms in 
the nanocontact were also multiplied by a random factor 
in the range 0-1.5. Averages were done over at least 
1000 realizations. To calculate the standard deviation, 
instead of directly working on the distribution dG/dV 
that gives very noisy results, we introduced the results for 
the transmission of each channel into Eq. . As shown 
in Ig both calculations agree in identifying the minima in 
fluctuations. The ECM results depicted in Fig [S] clearly 



show a peak in the conductance histogram near Qq, with 
a minimum in fluctuations. This means that only one 
channel is contributing to the unit conductance, being 
this value the most probable. The UHF results, however, 
show finite fluctuations over the whole range. 

Summarizing, the effects of electron-electron interac- 
tion on transport through hydogenated Pt nanocontacts 
have been investigated. Our results indicate that when 
two H are located transversally to the nanocontact axis 
[Sj, a Kondo effect stabilizes the conductance around 
one quantum. Spin-spin correlations are much like those 
which occur in the standard Kondo effect observed on 
two dots in a similar geometry [J| . The effect occurs due 
to the decoupling of d and s bands produced by the com- 
bination of the transversal location of the molecule and 
to the dxy symmetry of the Pt d-orbital. The result is 
robust and largely independent of the model parameters. 
Besides, a single channel contributes to the conductance, 
leading to a minimum in the conductance fluctuations, 
in agreement with experiments [l|. This analysis per- 
mits to conclude that the unit conductance observed in 
this system is the result of a combination of highly Kondo 
correlated spins and topology. 
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